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1 Introduction 



Combining the quantum theory and the general relativity is important for the unifi- 
cation of four elementary forces as well as for the smearing singularities appeared in the 
classical relativity and the exploring the fluctuation in the early universe. Black hole 
radiation is also one of important issues of quantum theory of gravity [|TJ. 

Canonical formalism of general relativity was formulated by Arnowitt-Deser-Misner 
(ADM) H and by Dirac ||. According to the invariance of general coordinate trans- 
formation, Hamiltonian dynamics leads to a totally constrained system. A canonical 
Hamiltonian treatment has been developed for Schwarzschild geometry by Berger et 
al. |3j] and Lund ||, by taking a special coordinate condition. They have shown that 
the reduced Hamiltonian describes no dynamics within their special choice of coordinate 
condition. This conclusion is compared to Birkhoff theorem in spherically symmetric ge- 
ometry ||. The Wheeler-DeWitt (WD) equation is obtained [0, by taking Schrodinger 
picture and by imposing constraints on the state vector. The Schrodinger picture is ad- 
equate to take into account of the space-time structure of the universe, though it does 
not incorporate with the unification of elementary forces, which is usually described in 
the Fock space formalism. 

There are serious problems of the quantum treatment, the definition of the time M, 



the interpretation of the wave function which is a solution of the WD equation|TO[, and 



the realization of classical system from quantum system[ll 



Since it is difficult to treat the WD equation in general, the spherically symmetric 
geometry is often considered. This geometry can describe the local black hole structure 
as well as the global cosmological structure, though it is not possible to include grav- 
itational waves. The Dirac quantization formalism and the WKB approximation has 



been applied to spherically symmetric gravity by Fischler, Morgan, and Polchinski [12 



They derived the WKB wave functions, and applied them to the nucleation of matter 
in the epoch of the inflation universe. 

The ADM formalism for the geometro-dynamics of the Schwarzschild black holes has 
been investigated extensively by Kuchaf |13[]. The reduced Hamiltonian was derived as 
the surface term imposing that the constraints held. The reduced Hamiltonian method 
has been developed for the Reissner-Nordstr6m-(anti-)de Sitter black holes by Louko et 



al., and Makela et al.p4|, fLo] . The Ashteckar's canonical formalism was also developed 
for the spherically symmetric gravity [[H], [T7|]. In Ref.JTB], an exact Hamilton- Jacobi 
solution of the (classical) WD equation was given. 

In solving the WD equation, there are operator ordering ambiguities and the regu- 
larization problem in defining the product of operators at the same space time point. 
DeWitt made the regularization scheme that the operator can commute at the same 
space-time point requiring 5(0) = 5(0)' = ||. Schwinger proposed a resolution of 
the operator ordering problem by using some symmetric ordering and taking the rela- 
tion 5(0)' = JIBfl . The operator ordering should be taken into account, but usually it 
is not taken seriously. If the ordering is neglected, the solution of the WD equations 
valid only in the semi-classical (or WKB) approximation. The relation between the 
operator ordering and the anomaly which arises from the algebra among constraints are 
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discussed by Kamenshchik and Lyakhovich in the analogy with string theory [Tj| . They 
introduced a new basis of gravity constraints and studied them by using the harmonic 
expansion for cosmological models. They got a relation between the space dimension 
and the spectrum of matter fields, just like the string theory. Recently Hori gave the 
solution of the Wheeler-DeWitt equation for the Schwarzschild-de Sitter geometry un- 
der a regularization that any operator commutes at the same space-time point, and 
considered the problem of time by using the Heisenberg equation of motion [^(J . 

In this paper, we consider the black hole solution in spherically symmetric geometry. 
We define the operator ordering for the Hamiltonian, the momentum and the mass 
operators such that they satisfy a closed algebra, and derive analytic solutions which 
satisfy all constraints. 

This paper is organized as follows. In section 2, we review and summarize the 
canonical formalism of the Einstein-Maxwell theory with a cosmological constant in 
spherically symmetric space-time. The derivation of the mass function is also described. 
In section 3, we perform the canonical quantization in the Schrodinger picture and 
obtain the set of constrained wave equations: the momentum constraint equation, the 
WD equation, and the mass and the charge eigen state equations. They are solved 
consistently after fixing the operator orderings. Discussion is given in section 4. 



2 Canonical formalism of the Einstein-Maxwell theory in spher- 
ically symmetric space-time 

In this section, we summarize the canonical formalism of the Einstein-Maxwell theory 
with a cosmological constant in the spherically symmetric space-time in four dimensions. 
We use the natural units c = h = G = 1 and follow mainly the convention adopted in 
Kuchaf's work[13]. 

We start to consider the general spherically symmetric space-time metric in the ADM 
decomposition 

ds 2 = -N 2 dt 2 + A 2 (dr + N r dt) 2 + R 2 dQ 2 , (2.1) 

where dfl 2 is the line element on the unit sphere, the lapse function N, the shift vec- 
tor N r , and metrics A and R, which are the function of time coordinate t and radial 
coordinate r. 

The action of the Einstein-Maxwell theory is of the form: 

I = -^Jd*x^g~(^R-2\-F 2 u ) , (2.2) 

where A denotes the cosmological term, and ^'R and are the scalar curvature and the 
determinant of metrics in four dimensions. The electromagnetic field strength is denoted 
by F^v, which is the function of only t and r. The only non vanishing component is 

F 01 = A 1 - A ', (2.3) 

where A and A\ denote the electromagnetic fields. Subscripts, dot and prime denote 
the derivative with respect to t and r, respectively. As the gravitational field and 
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the electromagnetic field depend only on the time and radial coordinates, this action 
describes the space-time structure of the universe, but does not include the gravitational 
waves and the electromagnetic waves. 

Inserting the metrics ( |2.1| ) into the action ( |2.2| ), the action in the ADM decomposition 
is given as 



■' ATr\2\ 



I = JdtJ dr[-N~ 1 (R(-A + (AN r )')(-R + R'N r ) + ~A(-i? + R'N 

+N(-A- 1 RR" - -A^R' 2 + A~ 2 RR'A') + -NA(1 - XR 2 ) 
2 2 

+ ^N- 1 A- 1 R 2 (A 1 -A' ) 2 ] . (2.4) 



The canonical momenta are obtained from the action fl2.4p by differentiating with respect 
to R, A, and A 1 : 

P A = —N~ 1 R(R — R'N r ) , 

p R = -JV _1 [i2(A- (AN r )') +A{R-R'N r )} , 

P A = N- 1 A- 1 R 2 {A 1 - A' ) . (2.5) 

The fact that the canonical momenta conjugate to N, N r and A vanish gives the pri- 
mary constraints. The action can be written in the canonical form after the Legendre 
transformation 

/ = JdtJ dr[P R R + P A A + P A A 1 - {NH + N r H r + A H A )} , (2.6) 

where 

H = -R~ 1 P A P R + ^R- 2 APl + A^RR" - A~ 2 RR'A' 

+ I A - 1 i?' 2 - I A + -AR 2 + -AR- 2 P 2 , 
2 2 2 2 A ' 

H r = R'P R -AP' K , 

H A = -P' A . (2.7) 

The quantities H and H r are called the super-Hamiltonian and the super-momentum of 
gravity. They should vanish 

H ^ , H r ^0 , H A ^0 . (2.8) 

The constraint for H A corresponds to the Gauss law of the electromagnetic field. 

The Poisson brackets at the same time t between these constraints form a involution 
relations 

{H r (r), H r (r')} = H r (r) 5'(r — r 1 ) — (r <-> r') , 
{H(r), H r (r')} = H'(r) 8{r - r') + H(r) 5'(r - r') , 
{#(r), H(r')} = A~ 2 (r)H r (r) 5'(r - r) - (r <-> r) , (2.9) 



3 



where r and r' are different space coordinate and (r <-> r') indicates the interchange the 
argument in the preceding term. The Poisson brackets between Ha and other constraints 
are trivially vanish. These constraints form the first class, that is, a kind of gauge theory, 
which reflects the general covariance of the spherically symmetric space-time. 

In order to study the canonical theory of black hole, we have to introduce the mass 
function. The mass function is defined by the integration with respect to the space 
coordinate r of a linear combination of the super-Hamiltonian and the super-momentum 

as 



= -U-R^Pl + A- 2 RR' 2 -R- R~ l P 2 A + \r*) - m 

= M-m, (2.10) 

where m is a integration constant, which is assigned a mass of black hole as we see 
later. The physical meaning of the mass function is studied by Nambu and Sasaki [2I| 



identifying with the energy of the inner part of the Schwarzschild black hole. The canon- 
ical treatment of the mass function was considered by Fischler, Morgan, and Polchinski 



12| , and by Kuchaf [13| for the Schwarzschild black holes. The mass function can be 
considered as a constraint, of which Poisson brackets with other constraints at the same 
time t are 

{M(r), H r (r')} = M\r)8{r - r') , 
{M(r),H(r')} = -A(r)- 3 R'(r)H r (r)5(r - r') , 
{M(r),M(r)} = 0. (2.11) 

The Poisson bracket shows that the mass function is a constant of motion in weak sense. 



3 Quantization, operator ordering and analytic solutions of 
Wheeler-DeWitt equation 

In this section we proceed the canonical quantization of the Einstein-Maxwell theory 
with the cosmological term in the spherically symmetric space-time. 

The momenta corresponding to the electromagnetic field A\ and the metrics A, R 
( |2.5| ) are represented by functional differential operators in the Schrodinger picture 

P A (r) = -i 



Pa(t) 



6A 1 {r) 
5 

5A{r) 
5 
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In the following we use notation "hat" for the quantized differential operators and we 
neglect to express the time t explicitly for operators because we treat them at the same 
time. The fact that momenta conjugate to Aq and N, N r and their time derivatives van- 
ish leads the constraint equations (|2.8|). In quantum theory, these constrained equations 
are understood as the equations operating on the state vector 

= , H r ^> = , H A m = 0, (3.2) 

and 

Mty = , (3.3) 

where the integration constant m plays a role of the eigen value of the mass operator. 
The first equations in Eq.( |3.2|) is famous Wheeler-DeWitt (WD) equation for spherically 
symmetric case. 

It is quite hard to obtain solutions of the above four constraints. Furthermore, 
we have to specify the operator ordering for H, H r and M to make the problem well 
defined. We utilize this operator ordering and choose it so that we can solve the above 
four constraint equations. Although we have a freedom to specify the operator ordering, 
it is a quite difficult problem to give a consistent ordering in a sense that H, H r and M 
form a closed algebra. 

Our strategy to obtain analytic solutions is as follows: 

(1) We solve the charge constraint Ha^ = and obtain general solutions. 

(2) We specify the operator ordering for H r and then obtain general solutions of the 

momentum constraint H r ^> = 0. 

(3) We define a specific operator ordering for M such that H r and M form a closed 

algebra. Then, we require the mass constraint (M — m)^f = for solutions of 
momentum constraint. 

(4) We define the Hamiltonian through the spacial derivative of the mass operator, 

M', which is defined to reduce to the classical one in Eq.(2.10) when the ordering 
is neglected. We derive an explicit form of H and show that solutions for three 
constraints Ha^ = 0, H r ^> = and (M — m)^ = satisfy automatically the 
Hamiltonian constraint H^> = 0. 

(5) We show that the operators H, H r and M form a closed algebra. 

(6) We take a specific operator ordering for M and then derive analytic solutions for 

all four constraints. 

Following the above strategy, we shall explain the procedure to solve the problem in 
detail. 

(a) General solutions for the charge constraint 
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The wave function \I/ is a functional of the electromagnetic field A\ and the metrics 
A, R , and is assumed to be in the form of separation of variables as 

= V A [A 1 ]V G [A,R] , (3.4) 

where ^a and ^g denote the wave function of the electromagnetic part and that of 
gravity part respectively. The Gauss law equation Ha^a = corresponds to the charge 
conservation and its eigen value equation is solved as 

Pa^a = Q^a , 

V A [A] = exp(* J dr QA x {r)) , (3.5) 

where Q is the eigen value of the conserved charge. 

(b) General solutions of the momentum constraint 
We defined the operator ordering of H r 



H r = R'P R - A(P A )' , (3.6) 
and the momentum constraint equation is explicitly written by 

Note that other form of operator ordering is related to that in Eq.(3.6) by the similarity 



transformation 22 



We define a quantity Z which is a functional of R and A and is independent of the 
coordinate r, and require that it satisfies the momentum constraint equation 

[Z,iT r (r)]=i(i?(r) I ^-A(r)(^)')=0. (3.8) 

The solution is expressed by 

Z = J drAf(R, A~ l R') = J dr J dAf(R, A _1 _R'), (3.9) 

where / and / is the arbitrary function of R and A _1 i?', and the integration of r extends 
between boundaries. The function / and / is related by 

r R'/\ rl T _ 

f(R, A-'R') = -A-'R* j —f(R, x) . (3.10) 

By using Z, the solution of the momentum constraint is given by 

^ G = ^ G (Z) } (3.11) 
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This solution is a general solution in a sense that the functional form of / or / and also 
is not specified yet. These forms are determined by imposing the mass constraint. 
We comment that the most general solution is expressed by the product as ^ G = 
■* G {Z U Z 2 ,...) where Z j = J drAf 3 (R, A^R'), (j = 1,2,...). 

(c) General solutions of the mass constraint 
We defined the operator ordering of M as 

M — m— ^R-'P^Pa - ^R( X - F) , (3.12) 



where 



P| A) = AP A A~\ 

X = A~ 2 R' 2 , F — 1 — 2mR~ 1 + Q 2 R~ 2 — —R 2 . (3.13) 

3 

The mass operator in Eq.(3.12) reduces to Eq.(2.10) when the operator ordering is 
ignored. The ordering factor A is taken as 

A = A z (Z)A(R,x). (3.14) 

The following comments are in order: For A, we assume that it is a function of Z, R 
and x — (A" 1 /?') 2 in Eq.(3.13). The reason of this choice is due to the fact that these 
operators as well as P\ are good operators in a sence that they satisfy the following 
commutation relation with H r 

[0(r), H r (r')} = i(j)'(r)5(r - r') , (3.15) 

where <fi = R, R' / A, P A and Z (Z' = 0) and a function of these quantities. This is 
needed for H r and M form a closed algebra as we shall see later. Then, we assume the 
factorization between a function of Z and a function of R and Xi which is needed to 
guarantee 

[A(r),Pf V)] = t5(r - r') , [pf ) '(r) , ) '(/)] = . (3.16) 

The mass operator M — m is applied to solutions in Eq.(3.11) of the momentum 
constraint. The mass constraint equation is now given by 



5Z\ 2 d 2 ^ G (Z) 
SA J dZ 2 



+ A 



5A\ 5A, 



dZ 



+ (RJ X - F) 2 ^ G (Z) = . (3.17) 



Now we fix the functional form of A defined in Eq.(3.14) and / = SZ/SA defined 
in Eq.(3.9) such that the constraint equation in Eq.(3.17) does not contain R and A 
explicitly, that is, it becomes an ordinary differential equation of Z. We take 
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Since / is fixed, the variable Z is determined. Now we only need to determine the 
functional form of ^a- The function \1/ G satisfies 

Only the unfixed part is A Z - We take the form of A z so that the equation can be solved 
analytically in below. 

(d) Hamiltonian 

The operator ordering of Hamiltonian is defined through the relation 

(M)' = -A^R'H - R^PPA- 1 ^ , (3.20) 
where this relation is a quantum version of Eq.(2.10) and 

P[ B) = \{Pa + Pl A) ) = A^P A A-^ . (3.21) 
The Hamiltonian is explicitly written by 

H = l -kR~ 2 P [ PP k - R^PrKP^K- 1 + AR'- 1 - F))' , (3.22) 

where 

PP = P[ A) - tRR'- 1 U- l5 ^\ = CP A C- 1 , (3.23) 



with 



C = A exp 



(3.24) 



In the limit of neglection opereator ordering, this Hamiltonian reduces to the classical 
one defined in Eq.(2.7). 

From Eq.(3.20), we see that solutions of H r ^G = and (M — m)^> G = satisfy 
automatically the Hamiltonian constraint H^g = 0. 

(e) Commutation relations 

For H r and M defined in Eqs.(3.6) and (3.12), we find 

[M(r), M(r')] = , [M(r), H r (r')} = iM'(r)5(r - r') , 



[H r (r), H r (r')} = i \H r (r)^6(r - r') - H r ^)-^S{r - r')j . (3.25) 

We comment that the mass operator M in (3.12) consists of good operators in a sense 
that we stated in Eq.(3.15), so that the commutator between H r and M satisfies the 
above simple and closed algebra as we expect form the Poisson brackets in Eq.(2.11). 
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The commutation relations involving H is complicated, but we can derive them by 
differenciating the above commutators among M and H r : 

= ^[M(r),M(r')} 

= [-A^R'H - R^P^A^Hr, M{r')\ . (3.26) 

We find 

[H{r),M{r')} = iA~ 3 R'H r (r)5(r — r') + A~ 2 R~ 1 H(r)5(0)5(r — r') 

1 A(r) 



4R'{r)R(r)R{r') 



-{g(r,r')P A (r') - P[ A \r')g(r,r')} A~\r)H r (r) , 

(3.27) 



where 



g(r,r') = — - \A 1 ( r ^^ A{ ' 



SA(r) \ y ' 5A(r') J 
- (A-/A' z )' mm + + l K *(0)} S(r - r>) 

= 9(r',r). (3.28) 

Here, the dash for Az denotes the derivative with respect to Z as A' z = dAz/dZ. As 
we expected, the main term iA~ 3 R'H r (r)5(r — r') coincides to the one appeared in the 
Poisson bracket (2.11). Other terms that are proportional to H and H r are those due to 
the operator orderings. Since the commutaror between H and M contains only terms 
proportional to H and H r , the mass operator M is the conserved quantity. 

By differentiation the commutator between M(r) and H r (r') with respect to r and 
by using the commutator between H r (r) and H r (r'), we find 

[H{r), H r (r')\ = iH'(r)5(r - r') + iH(r)5'(r - r') , (3.29) 

where 5'(r — r') = d5(r — r') / dr. The above commutator is just the quantum extensiton 
of the Poisson bracket in Eq.(2.9). 

By differentiation the commutator between M(r) and M(r') with respect to r and 
r' and by using the commutation relations obtained above, we find 

[H(r), H(r')} = iA~ 2 (r)H r (r)5'(r - r') - h(r, r')H r {r') - (r <-> r') , (3.30) 

where 



(3.31) 



The function h(r, r') is determined as a function of Az and A. We do not show the 
functional form of h(r, r') explicitly here, because it is complicated. The first term in 
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the right-hand side is what we expected from the expression of the Poisson bracket. The 
second term is due to our choice of the operator ordering. 

We showed all commutators and saw that the algebra is closed. The commutators 
correspond to the Poisson brackets in Eq.(2.9), except for the appearance of the terms 
which are proportional to H r . They contain g(r,r') and h(r,r') which are determined 
when A z and A are given, and are the non-local terms in the sense that they are not 
proportional to the delta function and its derivatives, unless we take A z l A' z to be a 
constant, i.e., A z = expaZ with a a constant. In below, we take A z some polynomial 
of Z, so that these non-local terms exist in the commutators. 

(f) Analytic solutions 

We choose A z such that the solution of the equation (3.19) becomes a special func- 
tion. 

(f-1) Hypergeometric type solutions 

If we choose A z = Z a (Z — I) 5 and transform 

^ { G ,5) {Z) = Z a+ \Z - l) s+1 ip(Z) , (3.32) 

we find 

„ f a + 2 S + 2 \ , . a + 5 + 2 , 

so that solutions are given by hypergeometric function as 

¥f\z) = Z^iZ-iy+^a.Fia^n-Z) 

+a 2 Z 1 - 7 F(a- 7 + l, / 3-7 + l,2- 7 ;Z)) , (3.34) 

where a\ and a 2 are integration constants, and a, (3 and 7 are given by 

a(3 = a + 5 + 2 , a + (3 = a + 5 + 3 , 7 = cr + 2. (3.35) 

(f-2) Bessel type solutions 

As a simpler case, we take A z = Z 2v ~ x . Then, the equation becomes 

d 2 ^ G {Z) 2v-l(M> G (Z) . , . 

Solutions of Eq.(3.32) are given by the Hankel (or Bessel) functions as 

*£>(Z) = Z» (h HW(Z) + b 2 Hi 2 \Z)) , (3.37) 

where b± , 62 are integration constants. In a special case of v — 1/2, solutions are simply 
expressed as 

^g /2) (Z) = J- (61 exp (i(Z - tt/2)) + b 2 exp (-i(Z - tt/2))) . (3.38) 
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4 Discussion 



We have studied the canonical quantum theory of the Einstein-Maxwell model with 
a cosmological term in spherically symmetric space-time. We defined the operator or- 
derings for H, H r and M and showed that they form a closed algebra. The commutators 
reproduced all terms appeared in the Poisson bracket. In addition, we found some terms 
which are proportional to H r . These terms contain in general the non-local terms in a 
sense that they are not proportional to 5 function and its derivatives. In this operator 
ordering, we derived analytic solutions which satisfy all constraints, the Hamiltonian, 
the momentum and the mass constraints. In our knowledge, this work is the first to 
obtain solutions of these constraints in the spherically symmetric geometry, with the 
operator ordering such that the algebra is closed. 

From the point of the quatnum field theory, the field operator in the the Schrodinger 
picture is generally difficult to define consistently and mathematically. For general rela- 
tivity, it is difficult to solve the Wheeler-DeWitt equation even formally. Our procedure 
succeeded in bypassing these difficulties through the mass operator, which is the integral 
of the super-Hamiltonian and the super-momentum. 

We now compare our solutions and the WKB (semi-classical) solution. The WKB 
wave function is given by (12], [16], |2(| 



~ ci exp (iZ WKB ) + c 2 exp (-iZ WKB ) 
where C\, c 2 are integration constants and 

J dr J A dA R^ X - F , 



J WKB 



J drRk(s] X -F + Vxkog | ij^rp I) > ( forf< °) 



where x an d F are defined in Eq.(3.13). The function Z WKB turns out to be the same 
as our Z in Eq.(3.9) and the WKB solution agrees with our Bessel type solution for 
v = 1/2. Our Bessel type solutions with arbitrary v approach to the WKB solution for 
large values of Z. 

Our solutions are rather general since we did not impose any coordinate conditions. 
All solutions that are known before are obtained by imposing some coordinate condi- 
tions. Thus, we limit the geometry such as Schwarzschild geometry or the de Sitter 
geometry, and impose some coordinate conditions, and see whether our solutions reduce 
to some known solutions. Since our solutions are given as a function of Z, this limiting 
procedure is the process to find the functional form of Z. 
(i) The Schwarzschild black hole case (Q = A = 0): 

Further, the coordinate conditions, R' = so that x = A~ 2 i?' 2 = 0, and R = 1 are 
imposed. Then, we require A = y/—F in addition, the function Z is determined by 

Z = J drkR^x ~ F | X=0 ,A =V ^F= / dr{-R + 2m) . 

The Bessel type solutions with the above Z agree with the ones given by Nakamura et 
al.[23] and by our previous work [24]. This case is extended to the case where Q, X ^ 0. 
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In our recent work[25], solutions are given and the interpretation of wave functions have 
been given. 

(ii) The de Sitter universe case (m — Q — 0): 

The coordinate conditions, x — 1 an d R = \J\/3R' are imposed further. If we 

require the relation A = ^J\/3R = a(t)e r , where a(t) = exp (y/~\/3t), Z is determined as 
Z oc a(t) 3 . If we choose the index of the Bessel type solution as v = 1/3 (Airy function), 
our solution agrees with that of Horiguchi[26]. 

In summary, our solution can describe the local structure such as a black hole and 
the global structure such as universe. 

Our solutions may be used to investigate the early universe in the inflation time. Also 
the solutions may be used for the analysis of the quantum universe and the quantum 
black hole. In order to study the quantum epoch beyond the semi-classical region, we 
assume some interpretation of quantum universe. The de Broglie-Bohm (pilot wave, 
quantum potential) interpretation[27] lies in an unique position in a sense that this 
provide the notion of trajectory in the fully quantum region. The study of the quantum 
black hole is now under investigation by using the de Broglie-Bohm interpretation. 
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